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Abstract 

It has been argued by Jacobson, Liberati and Mattingly that synchrotron radia- 
tion from the Crab Nebula imposes a stringent constraint on any modification of the 
dispersion relations of the electron that might be induced by quantum gravity. We 
supplement their analysis by deriving the spectrum of synchrotron radiation from 
the coupling of an electrically- charged particle to an external magnetic fields in the 
presence of quantum-gravity effects of the general form {E/Mqg)°'- We find that 
the synchrotron constraint from the Crab Nebula practically excludes a ^ 1.74 for 
Mqg ~ mp = 1.2 X 10^^ GeV. On the other hand, this analysis does not constrain 
any modification of the dispersion relation of the photon that might be induced by 
quantum gravity. We point out that such quantum-gravity effects need not obey the 
equivalence principle, a point exemplified by the Liouville-string D-particle model 
of space-time foam. This model suggests a linear modification of the dispersion re- 
lation for the photon, but not for the electron, and hence is compatible with known 
constraints from the Crab Nebula and elsewhere. 



CERN-TH/2003-197 
August 2003 



1 Introduction 



Modified dispersion relations have been first suggested in the context of stringy 
quantum gravity (QG) in 1 , based on the Liouville string approach to quantum 
space time j^I- Later, analogous modifications have also been suggested in the con- 
texts of other models of QG, either phenomenological |2] or theoretical. Examples 
of the latter include loop gravity |3] and novel models in which the Planck scale is 
viewed as a length that remains invariant under modified non-linear Lorentz trans- 
formations [5J. Violations of Lorentz symmetry had been suggested earlier [H] as 
a way of avoiding the Greisen-Zatsepin-Kuzmin (GZK) cutoff for ultra-high-energy 
cosmic rays (UHECRs). 

Liouville-string models of space-time foam |2j motivate corrections to the usual 
relativistic dispersion relations that are of first order in the particle energies, cor- 
responding to a vacuum refractive index ^ 1 — [E/Mqg)°''- a = 1. These effects 
are associated generically with deviations from conformal invariance in the effec- 
tive theory of low-energy excitations interacting with singular or topologically non- 
trivial quantum-gravitational (QG) degrees of freedom, inaccessible to low-energy 
observers |2] . Models with quadratic dependences of the vacuum refractive index on 
energy: a = 2 have also been considered jTj. 

The phenomenology of such models has grown rapidly. Following the original 
suggestion |Hlini to place bounds on the effective QG scale by comparing the arrival 
times of photons with different energies from gamma-ray bursters (GRBs), in order 
to probe the refractive index that might be induced by QG, it has also been pointed 
out that electrically-charged fermionic probes, either in astrophysics [10] or in ter- 
restrial atomic and nuclear physics experiments El E] , can constrain severely 
phenomenological models of space-time foam with modified dispersion relations. 

The most severe of all the known constraints seems to be that associated with 
synchrotron radiation from the Crab Nebula, whose sensitivity exceeds the Planck 
scale by (at least) seven orders of magnitude ^2); in the case of a linear modification 
of the dispersion relation for the electron. For photons, on the other hand, the most 
precise probe of a possible QG-induced refractive index comes from an analysis of 
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arrival times of emissions from GRBs 0, which impose Mqg ^ 10^^ GeV. Analyses 
of light from active galactic nuclei (AGNs) and pulsars J3] also reveal no evidence 
of QG effects, and may have comparable sensitivities to Mqg for photons. 

In the first part of this paper, we revisit the analysis of ^2]; by considering 
in detail the propagation of a charged matter probe interacting with an external 
magnetic field in the case of a generic modified dispersion relation, not necessarily 
linearly suppressed by the Planck scale. Our analytical result includes explicitly 
corrections due to the calculable change in the synchrotron radius in the presence of 
such QG effects. We confirm that the synchrotron constraint provided by the Crab 
Nebula on the electron's dispersion relation is robust theoretically. Unless it is 
relaxed for astrophysical reasons, this type of constraint is so strong that it already 
excludes an a < 1.7 correction to the vacuum refractive index for electrons, and has 
the potential to be sensitive to a quadratic correction: a = 2 in the near future, 
when higher-precision data become available. 

However, this does not mean that all types of QG corrections with a < 2 are 
excluded. The escape route is for QG to violate the equivalence principle, so that 
a = 1 and r] < 1 for photons, whilst rj = 1 for electrons. Remarkably, this is 
exactly what happens in the Liouville model of space-time foam proposed in [TK] . 
in the modern framework of the D-brane approach to QG jTHj. The reason for this 
violation of the equivalence principle for different categories of particles is explained 
briefly at the end of this paper, and is described in more detail elsewhere |17j . 

The structure of the article is as follows: in the next Section we summarize the 
Crab Nebula analysis of jT2j. Then, in Section 3 we present an analytical derivation 
of the synchrotron radiation spectrum in the presence of modified dispersion rela- 
tions, including the change in the synchrotron radius that we show to be a small 
effect. We also show that a quadratic modification of the electron's dispersion rela- 
tion with Mqg = mp is marginally excluded by the Crab data. Finally, in Section 4, 
we discuss the ability of this constraint to exclude some models of QG, showing how 
the specific model of space-time foam in ^H] evades it by violating the equivalence 
principle in a characteristic way. 
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2 The Synchrotron Radiation Constraint from the 
Crab Nebula 



Let us briefly summarize the main points of ^21; wlio first proposed considering tlie 
constraints implied by synchrotron radiation from the Crab Nebula. We highlight 
some steps which were not made explicit in their analysis, motivating our deriva- 
tion in the next Section of an analytical mathematical description of synchrotron 
radiation in the presence of modified dispersion relations. 

Following [12], we assume the modified dispersion relations (in units of the speed 
of light in vacuo c, which we now set to unity) 



for photons and electrons, respectively, where uj and k are the photon frequency 
and wave number, and E and p are the electron energy and momentum, with mo 
the electron rest mass. Here we assume linear QG effects, characterized by param- 
eters and ^e, extracting the Planck mass scale mp = 1.22 x 10^^ GeV. We also 
assume here that energy and momentum are conserved in particle interactions ^. 
We expect that the parameters and are negative semi-definite, given that in 
the framework of Liouville-inspired quantum gravity [15] there are at most sublu- 
minal modifications, if any, and superluminal modifications would cause troubles 
with gravitational Cerenkov radiation |^. We note that the modifications to the 
dispersion relation proposed in jT^j arise from a non- Minkowski induced metric in 
target space, which leads to some formal differences from the approach of ^2], as 
we comment later. 

We now recapitulate the derivation of constraints from the Crab Nebula syn- 
chrotron radiation on phenomenological models of quantum gravity that incorporate 
(|21) for electrons. It was assumed in J2| that the usual description of synchrotron 
^This may be questioned ^l^j, but any plausible violation is unimportant for our purposes. 



u^ik) 
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radiation was applicable, and the following facts about spectrum of the Crab nebula 
were used: 



• Synchrotron emission is observed up to energies of about 0.5 GeV, where the 
inverse Compton hump begins to dominate the spectrum. 

• Photons with energies up to 50 TeV are observed. 

• Electrons with energies up to at least 50 TeV are required by energy conser- 
vation to produce the observed flux of 50 TeV photons by inverse Compton 
scattering. 

In the standard Lorentz-invariant theory, the 0.5 GeV synchrotron emission in the 
magnetic field of the Crab nebula (~ 0.3 mG) is generated by electrons of energy 
5 X 10'^ TeV. The inference of this energy assumes, however, Lorentz invariance, 
which is just what we want to test. Hence the authors of [12j adopted the more 
conservative lower value of 50 TeV, which is inferred using only energy conservation 
in the inverse Compton process. 

In standard electrodynamics, accelerated electrons in a magnetic field emit syn- 
chrotron radiation with a spectrum that cuts off sharply at a frequency Uc which can 
be calculated as follows. According to the standard theory, electrons in an external 
magnetic field H , follow helical orbits transverse to the direction of if, with the 
orbit radius R given by 

(3) 

where E is the total energy of the electron. The orbital frequency of the electron in 
this case is: 

^0 = § (4) 

where (5± = v± is the component of the velocity of the electron perpendicular to the 
direction of the magnetic field. It was argued in [21] that there could be modifications 
to R that might affect the results of [12]. We show later that there are indeed QG 
induced modifications, but that these do not change the central result of 
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We recall that, in the standard Lorentz-invariant (LI) theory of synchrotron 
radiation, an electron moving in a magnetic field H emits a discrete spectrum of 
frequencies, which are integer multiples of: 

where 6 is the angle between the velocity of the electron, v, and the direction of H. 
The spectrum of the emitted radiation has a maximum at a critical frequency 

r r 3 eH 1 , . 

where e is the electron charge. The superfix LI in (jH)) stresses that this formula is 
based on a LI approach, in which one calculates the electron trajectory in a given 
magnetic field H and the radiation produced by a given current, using the relativistic 
relation between energy and velocity. All of these assumptions could in principle be 
affected by violations of Lorentz symmetry. 

It was stated in |T2] that, even without assuming Lorentz invariance, the critical 
frequency is given by 

= 3 ^ 1 

4R6{E) c{uc)-v{E) ^ ' 

where S{E) is the opening angle for the forward-directed radiation pattern, and c{uJc) 
and v{E) are the group velocities of the radiation and electron respectively. The 
self-consistent solution of ((Zj) for uJc{E) determines the cutoff synchrotron frequency. 
There was no attempt to solve this equation in J2] ■ Instead, it was argued that one 
can replace ciuc) by c = 1, so that the self-consistent solution is simply equal to the 
right-hand side of (jZj). To motivate this result, one may note that the electron and 
photon speeds are very close to the low-energy speed of light in vacuo c = 1. 

The synchrotron radius R for any given energy is determined by the equation of 
motion of the electron. The authors of ^2] assumed that gauge invariance is pre- 
served, and used the usual minimal coupling. To find the electron equation of motion 
in a magnetic field H, they used the dispersion relation Q for the Hamiltonian, with 
the momentum replaced by 

p ^ p - eA (8) 
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where A is a vector potential for the magnetic field. This yields the equation of 
motion a = [1 + ?)S,eE / 2M]{e / E) v x H, where only the term of lowest order in 
is kept, when one assumes E ^ mo ^. Since E <^ M, the authors of 12J argued 
that Lorentz violation would make very little difference to the orbital equation, 
so that the radius R would be related to the magnetic field H and the energy E 
of the electron by the standard formula Q- These arguments have been refuted 
qualitatively in [21J, with a subsequent reply in [22] • In order to arrive at a decisive 
conclusion, we embark in this article on a detailed quantitative analysis of the QG 
modification of the curvature radius. 

In the next Section, we revisit the analysis of T^, presenting a derivation that 
includes the modification of this formula for the orbit radius due to the modification 
of the dispersion relation. The potential importance of this modification has been 
stressed in |2I], but we find it to be relatively small. Our calculations should be 
accurate for all energies E ^ 0.3 Mp/\^e\- The analysis of |12j can then be used 
to infer a strong lower bound of the effective quantum gravity scale for electrons: 
Mp/\^e\ ~ 10^^ GeV. However, we reiterate that there are no such strict bounds 
for photons, for which the most stringent bounds on the quantum gravity scale 
Mp/I^^l ^ 10^^ GeV have been derived from gamma-ray bursts [0]. 

One might be tempted to assume universality for the QG corrections and 
^e, as a consequence of the equivalence principle, in which case the electron bound 
would exclude observable linear QG modifications of dispersion relations in general. 
However, one should not exclude a priori the possibility that the equivalence prin- 
ciple might be violated for such QG corrections, with different effective quantum 
scales for different particle species, as considered in [HI 123 ■ We argue in Section 5 
that this is indeed the case in the Liouville model for space-time foam of ^15ij, where 
electric charge conservation makes the foam medium transparent to charged probes 
such as electrons, but not to photons. 

Before proceeding with our detailed derivation of the modified synchrotron spec- 

^We note that the covariantization was done assuming that the background space-time is 
flat Minkowski space, as in the phenomenological analysis of jl2| . 
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trum in the next Section, we first recall how (3) was used in ^2] to derive synchrotron 
radiation bounds. Using the dispersion relations pi2p . and keeping terms suppressed 
only linearly by the QG scale Mp, one may write the difference of group velocities 
in the denominator of the last term of ((Tj) as: 

c{u)-v{E) = i^uo + {ml/2E^)-i,E, (9) 

The synchrotron radiation constraint comes from maximizing the electron group ve- 
locity with respect to the energy which yields for the maximum photon frequency 



Comparing (jlUp with ©, it was observed in [T2] that the factor mo'-f{E)/ E, which 
is different from unity as a consequence of Lorentz violation, is a bounded function 
of E, 

7(i^) = (i-.r^/^- (11-2^4)"''. (11) 

Maximizing Uc with respect to the energy E yields 

pH 

= 0.34 — (-eemo/M)-2/3, (12) 
mo 

which is attained at the energy E^ax = (-2mgM/5^e)^/^ = 10 {-^'^^^ TeV. The 
frequency tu™'^^ is the highest possible value of the critical synchrotron frequency 
for any electron energy. The rapid decay of synchrotron emission at frequencies 
larger than Uc implies that most of the flux at a given frequency in the synchrotron 
radiation peak observed coming from the Crab Nebula is due to electrons for which 
Uc is above that frequency. Thus cj™^^ must be greater than the maximum observed 
synchrotron emission frequency hhJohs-, which yields the constraint 



> . (13) 

mo V "^ot^obs / 



This was the main result of [12]. It implies, when one takes into account indicative 
physical values/estimates of the various quantities entering (fT^ . a very small upper 
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bound on the magnitude of the coefficient |^e| of the QG modification of the electron 
dispersion relation (which is assumed to be negative, corresponding to subluminal 
propagation) : 

l^el < 7 X 10-«. (14) 

As already mentioned, this analysis does not yield a stringent constraint on the 
magnitude of the photon coefficient |^^|. 

At this stage, we remark that, in the framework of our Liouville model of space- 
time foam, we have derived QG modifications only for massless neutral particles such 
as the photon and photino, the latter being used as a model for the neutrino |24j . 
As discussed later in this paper, different categories of particles may have different 
QG modifications in their dispersion relations, and these derivations do not apply 
to electrons. 



3 Analytic Derivation of the Synchrotron Radia- 
tion Constraint 

We now present an analytical derivation of the QG modification to the synchrotron- 
radiation formula, for generic subluminal modifications of the electron dispersion 
relation of the form: 

E'ip) = ml+p'-^ , (15) 

where Ai{=Mp/\C,e\ for a = 1) sets the effective QG scale in the modifications 
of the dispersion relation for the electron, which we seek to bound by means of 
the synchrotron radiation. We use the modified dispersion relation (fTH|) for the 
Hamiltonian of a charged spinless particle probe, which we place in an external 
magnetic field. It is convenient to parametrize the QG effects by means 

of the refractive index in vacuo rj: 

„=l-(-^)°, a>l. (16) 
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At ultra-relativistic energies, one may replace p° by E°' in the dispersion relation, 
and thus use the following for the Hamiltonian: 

E'=p' + ml-p'(^-^^ . (17) 

Notice that, in writing the dispersion relation in the form p7p . we have ignored 
terms proportional to the mass mo in the QG modifications, expecting them to 
be subleading, which is a self-consistent approximation at ultra-relativistic electron 
energies, as we shall see. Moreover, for the range of energies relevant to our problem, 
we can also ignore terms of order (E/Ai)'^" or higher in the following. 

Before analyzing the consequences of (fT7j) . we first make some important remarks 
about the Liouville foam model of fHI, as opposed to generic models of phenomeno- 
logical modified dispersion relations in flat space-times, such were considered in [12] . 
One important difference of the Liouville model of foam is that the modifications in 
the dispersion relation of a string probe owe their existence in the appearance of a 
non-trivial induced target-space metric: 

Goo = -1 ; Gij = 6ij ■ Go, = QsApi/Ms ~ ^gsP^/Ms (18) 

where Mg/gs = is the effective QG scale, Mg is the string scale, gs a string 
coupling, and Api is the momentum transfer during the scattering of the string with 
the D-particle defect jTH] • A generic dispersion relation for a particle with mass mg 
in the metric background ()18|) is: 

P.PuG^" = -ml , (19) 

which, as can readily be seen, implies jTHj a modified dispersion relation of suhluminal 
form (fTTj) . in the case where the direction of recoil of the D-particle defect is opposite 
to that of the incident particle, in the frame where the (massive) defect was initially 
at rest, 

P-^^~|^<0. (20) 

We recall that the requirement of subluminality emerges from specific properties of 
the underlying string theory, namely the Born-Infeld form of the effective action 
describing the dynamics of the excitations of the recoiling D-particle |15j . 
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Physically, this opposite recoil corresponds to capture of the particle probe with 
subsequent decay of the defect, and emission of a particle with modified momen- 
tum ^. It is this feature that may lead to violations of the equivalence principle, 
if electrons and photons have different probabilities for this process. This indeed 
happens in simple models, as we discuss later. 

We also observe from ()19|) that the inverse of the metric (jlSp results in a higher- 
order 'renormalization' of the rest-mass term rriQ in the dispersion relation (fT7j) : 
rriQ —>■ wig(l + We ignore such higher-order effects here, which is a self- 

consistent approximation for the range of energies relevant to the problem. 

In what follows, we keep our analysis as generic as possible, so as to minimize the 
dependence on specific models of foam. To this end, we use the modified dispersion 
relation (|17|) to define the momentum operator in a way that takes properly into 
account the QG modification ^: p ^ p = pr], which becomes p = p (l — in the 
linear case, where the minus sign indicates subluminal propagation. This definition 
implies that we define an 'effective' momentum squared as the part one has to 
subtract from the square of the energy in order to obtain the rest mass squared ^. 

It is this p that we shall couple to the external electromagnetic field. A, en- 
suring gauge invariance by minimal substitution, which implies the electromagnetic 
covariantization 

p^p+^A. (21) 

c 

As we now show in detail, in the case of QG-induced modifications of the electron and 

^Ref. |21| considered synchrotron radiation in terms of an interaction + -jh — > + 7, 
where denotes a soft photon representing the magnetic field. In our approach, we consider 
the quantum theory of synchrotron radiation |25) . in which an electron in a classical magnetic 
field emits a quantum of synchrotron radiation by changing the state of its wave function. This 
emission depends on a three-particle vertex that has been discussed in the context of our model of 
quantum gravity in jl9|, and could in principle violate energy conservation. However, such effects 
are important only when all three of the quantum particles have large momenta, which is not the 
case here, since the emitted photon is relatively soft. 

*In the Liouville foam model ^S], this definition could be thought of as the gravitational 
covariant derivative associated with the non-Minkowski induced metric (|18|l . 

^This definition resembles that in general relativity, where one defines the effective potential in 
a curved Schwarzschild space-time as the part one has to subtract from the square of the energy 
in order to obtain the square of the radial kinetic energy term. 
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photon dispersion relations, the maximal synchrotron frequency (jH)) gets modified, 
because of the modification of the orbital frequency as well as the velocity factor j3 
of the electron: 

where we now concentrate on the longitudinal motion. 

The modification of the orbital frequency is connected with QG effects on the 
radius of the orbit, due to the back-reaction of the foamy QG medium on the prop- 
agation of the electron. Such effects are purely quantum in nature, and therefore 
require a quantum treatment of the motion of an electron of charge e = — cq in 
a magnetic field, which we shall now present. For our purposes of estimating the 
change in the radius of the orbit as a result of the QG modifications of the dispersion 
relation, we do not take into account the spin of the electron, but consider instead 
the simpler problem of a spinless particle in a magnetic field. We justify the validity 
of this simplification later in the article. 

We consider a magnetic field along the z axis, which implies a vector potential 
Ax = —\yH , Ay = ^xH , Az = 0, and use the QG-induced dispersion relation 
(I17|l for the Hamiltonian operator of the quantum mechanical problem that describes 
the coupling of a spinless particle with this potential. The associated Klein-Gordon 
equation [23] can then be written in the following form (here we state explicitly the 
dependences on the speed of light in vacuo c) : 

{E^ -c\-p7] + -Ay -my}^(f,t) = , p=-thV, (23) 
c 

where rj is defined in flT^ . Proceeding, as usual, to the stationary case via 
ip{f,t) = exp{—i^t)ilj{r), we have: 

e - kl + rfV' -rj^A-V- ^A'^ i,{r) = 0, (24) 

where k = E/ch and ko = iriQc/h. In what follows, we consider only positive values 
of the energy, i.e., A; > 0, in which case we may rewrite in terms of cylindrical 
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A;2 _ kl + ri^Vl + r,^^ + - 2«r/7|- - 7V ) ^, ^) = 0, (25) 



coordinates {r,ip,z), as follows: 

where 7 = cqH /2ch. 

Following we now assume that the z component of momentum and the 
angular momentum are conserved as a consequence of rotational symmetry, which 
is maintained in the analysis of ^. Changing variables to ipi = ip/r] and z = z/t], 
we make the ansatz: 

V' = p^)I7iII7j/('-)- (26) 

where = 27Tn-^/L, and i and 713 are the azimuthal and vertical quantum numbers 
(corresponding to the appropriate components of the orbital angular momentum), 
which may take on positive and negative values, including zero. The radial part of 
the wave-function satisfies: 

^ dp^ dp 2 4 4/3772 J ^ ^ 

where p = 'yr'^/rj'^, and A = [k"^ — k^ — fc|)/47. 

The asymptotic behaviour as p ^ oo of the radial function f{p) can be chosen 
such that 

/oo = /(oo) ~ e-i , /o = /(O) ~ . (28) 
This allows us to write 

/ = foofouip) = e-(^/2) p'/' uip), (29) 
where u{p) satisfies the equation 



pu"+{i+l-p)u'+{X-e-l + ^il-v') + ^^{l-^,))u = 0. (30) 



^In certain models of anisotropic QG foam, angular momentum is not necessarily conserved. In 
such a case, our analysis here should be modified. 
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We now write this equation in a form that can be compared with the confluent 
hypergeometric equation 



pu" + (k — p)u' — ( a + flip H — ^ ) M = 0, 

p 



(31) 



where we set k 
the sul 
obtain 



+ 1, a = -A + £ + |, fli 



and as = f (1 - Making 
the substitution u = gp~'^^'^eP/'^ and changing the variable p = ^^^^^^^^ in (pTTjl . we 



4 pV4ai + 1 



p2 



(7 = 0. 



Substituting /c = (| — a)/ ^/Aai + 1, — w 4 2 ' 4 
form of a standard Whittaker equation 

k 



\ — t; + ^ + a?, we reduce 



-- + - + 
. 4 p 



1 2n 
p2 



(7 = 0. 



The solution for positive energies is 



1 



g = p™+i/2g-p/2^^^(^ ^l_k;2m + l; p), 



(32) 
to the 

(33) 
(34) 



where iFi{A,B,p) is the confluent hypergeometric function, which grows exponen- 
tially as p ^ 00. Therefore, the asymptotic behavior of the solution of (j30j) for large 
p takes the form 

(35) 



where T{x) is the standard Gamma function. This implies that, in order to respect 
the flniteness of the probability of the wave-function as p — > cxd, one has to impose 
the vanishing of the factor l/r(m + ^ — k), i.e., the selection rule: 

1 



m + 



k = —s 



5 = 0,1,2,. 



(36) 



where s is the well-known radial quantum number. Substituting the value of A given 
above, c.f., (^Tj) . we obtain: 



1,2 _ l2 _ l2 



47 



ST] + 



(37) 
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From this expression, we see that our analysis is vahd only in the regime of energies 
for which yf > 1/2, which implies 

{E/M)'' < 1 - ^ ~ 0.293 . (38) 
V 2 

This is not unreasonable, given that the concept of the effective field theory breaks 
down at the scale Ai. Within this range of energies, one can justify self-consistently 
the approximation of ignoring terms of order (i?/A^)^" or higher. 

We now compare the QG-modified spectrum (jTTjl with that appearing in the 
standard Lorentz-invariant theory, namely: 

+ + (39) 

where h = 0, 1,2, . . . is the principal quantum number. In the relativistic case of 
interest to us, both the spectra (jH^ and (jTfjl become quasi-continuous, since the 
quantum number n becomes very large. QG corrections simply introduce a small 
rescaling (by a factor of order one) of the eigenvalues of the energy states which 
define the quantized orbiting radii. 

Assuming that the motion is in the plane of the orbit {k^ = 0), one can then 
determine the radius R of the orbit just by equating the energy eigenvalues as given 
by (pITjl with those obtained in the case of a particle in a magnetic field: 



ch 

from which we obtain 



1/2 

= eoHRQG, (40) 



1/2 



rqc- r ^ " ^ ^ — ^1 . (41) 



for the QG-modified radius Rqg- 

We are interested in macroscopic values of the orbit radius, which, in the Crab 
Nebula case discussed in jT2], lie in the range ~ 10^°cm. Comparing this with the 
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other scales in the problem, we see that we are in a situation where n >> 1. Thus, 
the equation (jlUj) can be approximated by 

' ' '2^1^''' , (42) 



Ro V2\ V V^y 
where 

is the radius in the Lorentz-invariant theory. 

Some comments about this expression are now in order. First, we remark on 
the functional dependence of Rq on the principal quantum number n. In standard 
quantum mechanics, the mean square radius for a particle of spin s in a magnetic 
field is given by: Rq = I '^fi,s,^'^n,s,^rd^r, where is the wave function, and 
the sum is over possible spin states ^. In the case of the Klein-Gordon equation, Rq 
is obtained by solving (j3Up exactly in the limit f] = 1, with the result 

/ ~ \ 1/2 / I 3 



On the other hand, in the case of a particle with spin s = 1/2, corresponding to the 
Dirac equation in an external field, one obtains 

/ ~ \ 1/2 / ,1 
-Dirac I ^ \ I ^ . ^ ' 2 



^» =1-1 (l + ^l- («) 



If one is interested in macroscopically large trajectories, as is the case for ^2] and 
ourselves, where n ^ 1, we see from (I44|l . and that all three formulae give 
approximately the same result. In particular, the spin dependence of the shape of 
the orbit is not significant. This justifies our initial simplification of working with a 
spinless particle instead of a Dirac fermion in order to estimate the QG modifications 
to the electron orbit. 

This analysis demonstrates that QG modifies the critical synchrotron frequency 
fl22|) for radiation along the forward direction, in two ways. First by changing the 
relation of P to energy, due to the QG modification of the dispersion relation, but 
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also by modifying the orbit of the particle in the magnetic field, by slightly shrinking 
the average radius (021) • Thus, we may parametrize the maximal frequency © as 
follows: 

^c^^' = "t: — ^^^, (46) 



QQ _ 3 eH 1 Rq 



where the dependence on Rqg is dictated by the assumption of angular momentum 
conservation - the smaller the radius the higher the angular frequency. It is impor- 
tant to notice that, once one accepts the dispersion relation (fTTj). the result pi)) 
should be considered as exact for the entire range of energies for which ()H8|1 is valid, 
up to terms of higher order in (E/Ai)'^. 

Using the modified dispersion relation, then, identifying P^'^ with the group 
velocity, and keeping only terms linear in (E/Ai)"), we obtain the following result 
for the group velocity: 



dE p fa \ f E Y 



l-/5Sc-ll + (« + l)(^)". (47) 

On account of ()42j) . then, this yields for the QG modification to the critical syn- 
chrotron radiation frequency: 

' 72 mo (1 + ^/23W)V2 {§ + ia + 1) (A)") ' 

This function is plotted schematically (for a = 1) in Fig. ^ 

We observe that, since the effects on the radius are small for all values of 
E / M. < 0.3, the discussion of |12j reviewed in the previous Section remains qual- 
itatively correct. In particular, there is a global maximum within the range of 
energies allowed by namely i^^^ax' which should be higher than the ob- 

served synchrotron energy of 0.5 GeV. The maximum occurs for electron energies 
Emax = {2mlM''/a{a + 1))^^^^+"^ From this, one may obtain bounds on a, if one 
sets = Mp ~ 10^^ GeV, or, alternatively, obtain bounds for |^e| if one sets a to 
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Figure 1: The scaled critical frequency for synchrotron radiation in quantum gravity: 
u^^/ujc (uJc = 3eH/2mo) as a function of the electron energy E in units of M., as 
obtained from for a = 1. The synchrotron frequency is a bounded function of 
E, as in lijjj /. a feature which is not affected by the QG modifications of the orbit. 



a fixed value. It can be easily seen that the general formula (up to terms of order 
one) for the upper bound on |^e| is: 



'^'^'^WoJ UoJU(«+i)J U+2J ■ ^^^^ 

If the resulting bound on |^e| is less than unity, the sensitivity exceeds the Planck 

scale for the given value of a. For the linear case a = 1, the existence of the 

maximum implies that the discussion of ^2] remains intact and one arrives at the 

bounds of given in (fT3j) . (fT^ . 

It is easy to see that for the average magnetic field of Crab Nebula ^ Hcons = 

^It should be emphasized that the estimate of the end-point energy of the Crab synchrotron 
spectrum and of the magnetic field used above are indirect values based on the predictions of 
the Synchrotron Self-Compton (SSC) model of very- high-energy emission from Crab Nebula |27| . 
We use a choice of parameters which gives good agreement between the experimental data on 
high-energy emission and the predictions of the SSC model |57| EH ■ 
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260/iG one obtains < 6.8 ■ lO^^ (8.8 ■ IQ-^i)"- , which imphes that a < 1.72 
is excluded. If the lower value for Hncons = 160/iG |2Ij is used instead, then |^e| < 1 
for a < 1.74; we also observe that for a = 2 |^e| < (9(30 — 60) for the range of the 
magnetic field considered above. These imply already a sensitivity to quadratic QG 
corrections. Therefore, 

a > etc : 1.72 < < 1-74 (50) 

is the phenomeno logically allowed range of a, where the lower (upper) limit corre- 
sponds to Hcons (Hncons) ■ 

We should note that there are still some ambiguities regarding the magnitude of 
the maximal energy of synchrotron radiation before transition to inverse Compton 
emission, which could range from 30 MeV to 0.5 GeV. This uncertainty should be 
resolved by future gamma-ray detectors such as GLAST 13j, which could provide 
a better determination of the unpulsed gamma-ray spectrum in the energy range 
above 30 MeV, and thus a more precise determination of the maximum electron 
energy in the Crab Nebula. 

Measurements of high-energy emissions from the Vela pulsar jSHI indicate the 
operation of a mechanism for the production of very-high-energy gamma rays similar 
to that in Crab Nebula. This leads to an estimate of the magnetic field of order 
Hveia ~ 3fiG. When used in (Hnjl . this yields ac = 2.04, under the assumption that 
the first synchrotron hump is at the same location as in the Crab Nebula, namely 
0.5 GeV. If, one the other hand, the hump is at 30 MeV, then ac = 1.86. Therefore, 
improvements in measurements of synchrotron emission from the Vela Nebula (or 
other similar sources), which may be achieved in the not-too-distant future, could 
bring the sensitivity of such experiments closer to testing quadratic modifications of 
the electron dispersion relation. 

This will be decisive for probing models with quadratic suppression of QG ef- 
fects in the dispersion relation for electrons. However, we remind the reader that 

^If the end-point of the Crab synchrotron spectrum is as low as 30 MeV, then the upper and 
lower limits for ac in H50|) become 1.56 and 1.58 respectively. 
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Figure 2: In the model of Liouville foam of [IS], only string particles (S) neutral 
under the (unbroken) standard model group can he captured by the D-particle defects 
(D) in space time. This results in modified (suhluminal) dispersion relations for S, 
as a consequence of the recoil of the defect D in a direction opposite to the incident 
beam, after the capture stage. 

synchrotron radiation provides no significant bounds for QG modifications to tlie 
dispersion relation for photons [,12| . 

4 Violations of the Equivalence Principle in the 
Liouville Model for Space-Time Foam 

As already mentioned, the synchrotron radiation constraint 12 on the dispersion 
relation for electrons, refined in the previous Section, does not exclude the possibility 
of a larger QG modification of the dispersion relation for photons. Of course, this 
would require the violation of the equivalence principle for energetic particles. Re- 
markably, this is exactly what is predicted by our Liouville string/D-particle model 
for space-time foam, according to which only gauge bosons such as photons might 
have QG-modified dispersion relations, and not charged matter particles such as 
electrons. 

This difference may be traced to a cornerstone of D-brane physics, namely that 
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excitations which are charged under the gauge group are represented by open strings 
with their ends attached to the D-brane and only nenira/ excitations are allowed 
to propagate in the bulk space transverse to the brane. 

Our D-particle model of QG foam is based on point-like defects in space-time, 
that are nothing but zero-space-dimensional D-branes, embedded in a four- (or 
higher-) dimensional bulk Minkowski space-time. These affect the propagation of 
closed-string states via a recoil process [T3] (see Fig. that is associated with 
back-reaction effects in target space, which lead to the QG-modified effective metric 
(US)). In this picture, a closed-string state propagating in the bulk can be captured 
by the defect and split into two open-string excitations with their ends attached 
to the defect, which later recombine to become a closed-string state. According 
to the above-mentioned property of D-branes, the open-string excitations of the 
D-particle defect result in massless U(l) excitations for a single defect, or massless 
gauge bosons in the adjoint representation of U(N) for a group of coincident recoiling 
D-particles [SHI EH- 

In the original formulation of this D-particle model of QG foam pH], the D- 
particles have no 'hair', i.e., they have only vacuum quantum numbers. Therefore, 
they can absorb closed-string states only if they also carry no conserved charges. 
This would exclude any QG medium effects on particles with electromagnetic or 
colour charge, certainly including electrons and quarks. On the other hand, a QG 
effect would be expected for the photon, which carries no conserved charge. Because 
the electron has no interaction with the QG vacuum medium in this approach, it 
emits no Cerenkov radiation, despite traveling faster than photons, thus avoiding 
the vacuum Cerenkov radiation constraint considered by j32| . 

Neutrinos are mainly doublets of the SU(2) subgroup of the Standard Model. 
However, in generic seesaw models of neutrino masses, the light neutrinos mix with 
0{mw/M) components of gauge-singlet states, perhaps opening the window for 
QG effects, depending on details of the model. However, any such effects would 
be suppressed by at least a similar factor. In the case of gluons, which carry non- 
Abelian charges, it is possible that ensembles of D-particle defects would have a QG 
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effect on their dispersion relations, which would in general be unrelated to that on 
the photon. 

However, such a relation might appear if the U(l) of electromagnetism is actu- 
ally embedded in some simple GUT group such as SU(5) j^Sl, so that the photon 
becomes some combination of non-Abelian gauge fields, in which case QG effects on 
its dispersion relation might be related to those of the gluons. On the other hand, 
there exist non-simple GUT models, such as flipped SU(5) x U(l) jSl], in which the 
photon contains a significant component from outside the non-Abelian group factor, 
in which case the photon's dispersion relation will in general differ from those of the 
gluons. 

Even if the dispersion relations of individual gluons are modified by QG effects, 
these may be suppressed for colour-singlet hadrons such as protons. A QG modifi- 
cation of the proton's dispersion relation has been suggested as a way of explaining 
the possible existence of UHECRs beyond the GZK cut-off ^U], where energy non- 
conservation in high-energy reactions might also become significant fH^ ITH] . As 
observed in in the case of a linear QG modification of the proton's dispersion 
relation similar to (2), the magnitude of the analogous parameter would need to 
be of order lO^^^ - lO^^^ Such a suppression is (coincidentally?) 0{Aqcd/Mp), 
and such a suppression cannot be excluded at present, though a lot more work is re- 
quired before such a conclusion could be reached in the framework of the D-particle 
model of However, the above discussion makes it clear that the details of the 
dynamics of the interactions between particles and space-time foam are potentially 
important, and any naive phenomenological assumption of the equivalence principle 
could be misleading. 

Before closing, we comment briefly on recent claims that modified dispersion 
relations for photons would result in phase incoherence of light from distant galax- 
ies ISi]. Based on this suggestion, these authors proposed a stellar interferometry 
technique, and claimed to place stringent bounds on the effective quantum-gravity 
scale for photons, excluding linearly-modified dispersions. However, their arguments 
have been criticized in ^36,, on the basis that they overestimated the induced inco- 
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herent effects by a large factor (L/A)", where L is the distance of the source, A is 
the photon wavelength, and a is the parameter in the modified dispersion relation 
of the photon probe. The correct amount of cumulative phase incoherence induced 
by quantum gravity is jSH]: ^'P ~ (-^Z-^)^ Thus, for linear dispersion relations for 
photons (a = 1) such as those proposed in [Hllini, the technique of ^| cannot be 
used. Moreover, in the case of the D-particle recoil model for Liouville foam the 
re-emission of the photon from the recoiling D-particle (c.f. Fig. ^ is accompanied 
by a (random) phase in the photon's wave-function. This would destroy any cumula- 
tive phase incoherence, in contrast to the claims in Thus, a linear modification 
of the photon dispersion relation cannot be excluded by this argument. 

Our analysis emphasizes the interest in probing independently the dispersion 
relation of the photon. The study of the arrival times of photons from gamma-ray 
bursts still appears to be the best experimental probe of any possible refractive 
index for photons, and should be pursued further in the future. 

The mechanism for violating the equivalence principle discussed in this article 
may not be the only route to consistency with the current experimental constraints. 
However, it is certainly one promising way, in the sense that it does not invoke 
exceedingly small parameters, and appears naturally within a class of stringy mod- 
els of quantum gravity [E]. In other approaches to quantum gravity, such as loop 
canonical quantum gravity [Ij, one may encounter models which feature linear mod- 
ifications in the dispersion relations for matter excitations that are characterized by 
extra small parameters which can be bound by experiments. Such a scenario has 
been discussed recently in jSZ], where the structure of canonical (quantum) commu- 
tators between momenta and position operators in loop-gravity models is modified 
by extra scaling terms, described by a set of coefficients that are essentially free pa- 
rameters of the model. Such extra terms lead to linear modifications in dispersion 
relations, but with coefficients that are proportional to these parameters. The latter 
can then be bounded by requiring consistency of the models with the current exper- 
imental situation. In our opinion, the disadvantage of such an approach lies in the 
unnaturally small values one obtains for such parameters in this case, which is to be 
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contrasted with our case above, where the numerical coefficient in front of the hnear 
modification of the photon dispersion relation is naturally of order one. However, 
we cannot exclude the possibility this, or some other way yet to be invented, may 
be the way that Nature evades the current experimental constraints on quantum 
gravity. 
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